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TUn SUf-FUR BUOGBT OF THR TROPOSPHERG 


By 

S.N. Tiwari^ ami T.R. Augustsson^ 


SUMMARY 

A one-dimensional photochemical ti'opospheric model has been used to 
calculate the vortical profiles of tropospheric species. Particular 
attention is focused on the recent inclusion of the chemistry of the 
sulfur group, which consists of 15 species involving a total of 45 
chemical reactions. It is found that the chemistry of the sulfur species, 
because it is laa’gely anthropogenic, plays an increasingly important I'ole 
in the distribution of tropospheric gases, The calculated vertical 
profiles are compared to available measurements and generally found to 
be in good agreement. 


I.MTRODUCTION 

During the last few years there has been a growing concern over 
possible anthropogenic effects on the troposphere. The chemistry of 
the sulfur group is particularly interesting since the omissions of the 
species in the sulfur family are to a large extent anthropogenic in origin. 
Of all the tropospheric chemical families, sulfur is perhaps the least 
studied and understood. Thus, a detailed scheme to calculate the vertical 
distributions of 13 sulfur gases in the troposphere has beer, deveioped. 
Studies of sulfur in the atmosphere have primarily been concerned with the 
so-called "aerosol layer" (refs. 1 - 5) or concentrated n highly specific 
reaction paths and mechanisms in the troposphere (refs. 0-7) or certain 
regions (refs. 8-10), Very few global assessments of sulfur species 
exist (refs. 11 - 12), and those studies are qualitative in nature. 

^ Eminent Professor, Department of Mechanical Engineering and Mechanics, 

Old Dominion University, Norfolk, Virginia 23508. 

^ Graduate Research Assistant, Department of Mechanical Engineering and 
Mechanics, Old Dominion University, Norfolk, Virginia 23508. 



The recently developed sulfur package allows for quantification (ref. 13). 

The research work in this area has been going on for some time, and 
we now have a one-dimensional, time-independent, global, tropospheric, 
photochemical model. The model curx’ently includes the chemistry and 
photochemistry of OH, HO^, NO^^C* NO + NOg) , HNO^, NO3, NoOg, HNO2, CH3, CH3O 
CII3OOH, CH3O, H2CO, HCO, 0, 0('D), O3, H, .Nil, .NHa, NII3, N2H2, N2H3, 

,N2ll4, S, SO, SO2, SO3, HSO3, H2SO4, HS, H2S, CS, CS2, CH3SH. CH3SCII3, 
and COS. Of those 36 species, 16 are being transported (i.o., wo use 
the continuity equntio.n) and the remairung 20 are assumed to be in photo- 
chemical equilibrium (because they are short lived). There are 108 
reactions and 12 photodissociations in the cui’rent model, and the vertical 
transport is simulated by the use of eddy diffusivities. The results ot 
the present model have been compared with i*esults of other models with 
similar chemistry, x’ate coefficients, and eddy diffusion pax'.nmeters . 

In general, good agreements between various results has been found. 

Recently, an "ammonia package" has been included in the present model. 

This consists of the chemistiy of NH3 and five other closely related 
species. A total of 21 chemical reactions has been added to describe 
accurately the vertical profiles of the "ammonia group." Ammonia has 
been the subject of intensive studies, and the results luave been reported 
in references 14 to 17. 

The significance of coupled radiative-chemical studies of the Earth's 
troposphere has been emphasized strongly in the NASA Reference Publication 
1022 titled "Man's Impact on the Troposphere" ( Lectures in Tropospheric 
Chemistry , ref. 11). It is evident from this publication that the chemical 
composition of the troposphere is determined by complex geochemical 
cycles of carbon, nitrogen, sulfur, oxygen, and hydrogen species between 
the atmosphere, the oceans, the crust, and the biosphere. In recent 
years, anthropogenic activities have become a significant source of tropo- 
spheric species. As stated in the above reference, changes in the natural 
abundance of tropospheric species may alter other species via chemistry 
as well as produce perturbations in climate. For example, tropospheric 
ozone is strongly affected by tropospheric levels of nitrogen dioxide 
and methane, and stratospheric ozone is controlled by tropospheric nitrous 
oxide. Many of these anthropogenically produced species have absoi'ption 


lines within the 8- to 12-yin atmospheric window, which contributes to 
the greenhouse effect. These greenhouse absorbers include carbon dioxide, 
methane, niti’ous oxide, amjnonia, and ozone. The information available 
in the cited reference, therefore, indicates the need for coupled I'adiatlve- 
chemical studios of the Barth's troposphere in a systematic manner. 

The coupling between the Naval Research Laboratory's radiative 
transfer model (which includes Rayleigh and Mio scattering and surface 
albedo) to the NASA/LaRC one-dimensional tropospheric photochemical model 
has recently been completed. The preliminary results indicate that the 
vertical distributions of several key tropospheric species are altered 
significantly (ref. 18). The next report will address this topic in 
more detail. 

The next section of this report ("Tropospheric Sulfur Species") 
reviews and summarizes the existing literature on sulfur species. Each 
sulfur species is discussed briefly and available measurements are 
identified. The section titled "One-Dimensional Tropospheric Photochemical 
Model" then discusses the philosophy of modeling as well as the NASA/LaRC 
one-dimensional global, tropospheric, photochemical model used for the 
calculations summarized under "Results and Discussion." This model now 
consists of 37 chemical species (13 of which are in the sulfur group) 
and 108 reactions (of which 45 are sulfur related). The major text of 
"results and Discussion" describes all the reactions, productions, losses 
and lifetimes for each sulfur species as calculated by the present work. 

We have found that the sulfur species family consists of gases in various 
states of oxidation (from highly reduced to highly oxidized) with a wide 
range of lifetimes (from 10"'^ s to 4 mo) . In the same section we identify 
the various sources (i.e., natural vs. anthropogenic) and sinks (homogeneous 
vs. heterogeneous) for all sulfur species. In the concluding section 
("Future Perturbations to the Sulfur Budget"), some global implications 
resulting from perturbations in the sulfur budget are discussed. The 
continuity equation and its numerical solution together with appropriate 
boundary conditions are discussed in Appendixes A to C. Appendix D 
describes the calculations of incident solar radiation and photodissocia- 
tion rates, while Appendix E discusses the convergence criteria employed 
in this study. 


TROPOSPHERIC SULFUR SPECIES 


Introduction 


According to R. Stewart ot al. U» P- ddl : 

The abundances, sources and sinks, and photocl»omistrv of 
atmospheric sulfur compounds are all poorlv understood relative 
to oxygen, carbon or nitrogen compounds. Various sulfur species 
are important components of air pollution in many parts of the 
world. The anthropogenic input to the sulfur cycle is clearly 
significant and likely to increase in magnitude, but the con- 
sequences of these facts have been relatively little explored. 

They further conclude (p, 4b) that: 

Man is also intervening in the sulfur cycle, but global- 
scale consequences of this influence have yet to be identified. 

Host of the sulfur dioxide that enters the atmosphere is of 
industrial origin end can potentially alter such things as the 
atmospheric aerosol loading and the acidity of rainfall over 
at least mesoscale regions. The 26 Mt/yr of sulfur added to the 
soil in fertilisers is relativelv large compared with natural 
sources. If substantially increased coal burning occurs in the 
future in response to shortages of alternate fuels, the magnitude 
of the industrial sulfur source will increase still further, It 
is essential that some attempts to understand the potential impact 
of this fact be made in future research. 

The purpose of this study is to e.xplore and try to understand the 
implications that an increase in sulfur loading will have on global ti'opospheric 
chemistry. It has been known for a long time that species of the sulfur 
family adversely affect the air quality on the local and regional scales. 

This is due mainly to anthropogenic emissions in lai’ge metropolitan areas, 
'furthermore, it has been known for about two decades that species of the 
sulfur family diffuse upward into the stratosphere whore they form sulfate 
particulates tliat constitute the so-called "aerosol layer" usually found at 
around 20 km (ref. 21. Natui'al omissions of particulates from volcanoes 
also contribute to the aerosol layer. 

The sulfur gx'oup consists of species in various reduced and oxidised 
states and with widely varying lifetimes. All sulfur species have in common 
an unrelenting trend towax'd further oxidation. The oxidation chain for the 
sulfur family is shown schematically in figure 1. 
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Sulfur Dioxide CSOi2 “* mixing ratio o 5.0 x 10“ 

In largo concentrations, SO 2 is x^ty toxic to animals and plants. 

It is an important precursor to formation of p.rffjieuiates. This affects 
the aerosol loading, which has climatic implications in addition to the very 
obvious problem of reduced visibility, Another climatic aspect to consider 
is the thermal infrared absorption bands that SOg has at 8.? nm and 7.S nm 
iref. 19). At the present time the measurements of SOg are very inconclusive 
and in some cases contradictory. For example, Georgii and Uravenhorst 
tref. 20} and Maroulis et al. tref. 21) found that sulfur dioxide increase.s 
with altitude, while Georgii (I’ef. 22) found a decrease with altitude. 

This fact is currently a source of some controversy in the scientific 
literature. 


Carbonyl Sulfide (COS — mixing I'atlo » 5.0 x 10“^^) 

COS is a very stable molecule and is consequently well mi,xcd in the 
troposphere. It dl^^nisos into the stratosphere where it is believed to 
undergo photodissociation. It subsequently recombines with other species 
to form SOg, which, in the stratosphere, is believed to be the gaseous 
precursor of the aerosol layer (refs. 1-2). The observational data of 
COS indicate a relatively uniform distribution. Sandal Is and Penkett 
iref. 23) found values varying from 4.0 x 10"*'^ to S.O x 10“^^ with a moan 
concentration of 5.1 x 10 Maroulis et al. (ref. 24) measured averaged 
concentrations of 4.35 x lO"^® in Philadelphia, 4. 54 x lO'iO at Wallops 
Island (Virginia), and 5.11 x 10“^® at Lawton (Oklahoma). Torres et al. 

(ref. 25) found values ranging from 4 to 5 x 10"^^ during the so-called 
"GAMBTAG" flights. The average value during those flights was 4.7 .x 10"^*^. 

Hanst et al. (ref. 26) found values of about 2 to 3 x in some measurements 

that had largo e.xpei'imontal uncertainties associated with them. 

Hydrogen Sulfide (H^S — - mixing ratio = 2.0 x 10”^®) 

Hydrogen sulfide is formed primarily by microbial activity in the soil 
and in tidal flats (ref. 27), but it is also emitted by volcanoes (ref. 28) 
as well as by some industrial activities such as wood pulping (ref. 29) and sewage 
treatment (ref. 30). Measurements of the vortical distribution are very 
sparse, but seem to indicate a very rapid decrease with altitude, particularly 
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in the first few meteiT. Crof. 22). The mi.'dng ratio at the surface has 
been estimated at 0,2 parts per billion by volume (ppby)(rof, 31). The 
measurement*’ by Georgii (ref. 22} seem to verify this value. 

Carbon Disulfide (L'Sg mi.xing ratio « 1.0 x 10"^®) 

CS2 is formed both naturally by volcanic emissions {ref. 32) and by 
anthropogenic emissions such as manufacturing of petroleum fref. 33) and 
of synthetic fiber (E.P.A,, 1973-107S, refs. 34-36), CS« has been measured 
in seawater (ref. 37), giving an almost constant value of 5.2 x 10"^^ 
g/ml of water over a wide range of latitudes. This would indicate a natural 
source for carbon disulfide. Atmosphoi’ic measurements vary from 7 x 10”^^ 
to 3.7 X 10"^® with an average value of 1.9 x 10“-*5 (ref. 23). It has 
recently been suggested in references 4 and 6 that the oxidation of CS,7 
to form SO2 might bo similar to the oxidation for COS, It has also been 
argued that CS^ might oxidise to form COS (ref. 4), In either case, CSn 
is on important precursor for aerosol formation. 

Sulfuric Acid (HgSOi, — mixing ration 1,0 x 10“^0) 

H2SO4 is formed chemically as well as naturally, by volcanic omission 
(ref. S8), and anthropogenically, in automobiles a*ef. 36) and in H2SO4 
manufacturing (ref. 39), Sulfuric acid occurs both in liquid and solid 
forms. It forms the solid form (sulfate, SO^“) when it recombines with 
cations other than hydrogen. Hence, H2SO4 affects the environment in 
two ways; 

(i) by formation of aerosol particulates, which reduce both incoming 
radiation and visibility, and 

(ii) by lowering the pH in rainwater. 

In Scandinavia, which is downwind from the large industrial areas in 
Germany and Great Britain, a gradual lowering of pH in rainwater has been 
observed. The average rate of decrease has been 0,3 to 0.4 units per decade 
(ref. 40). Similar acidic rainfalls have been observed in the Atlantic 
provinces of Canada, downwind of the highly industrialised metropolitan 
"northeast corridor" (ref. 9) . Acid precipitation is kno\m to have an 
adverse effect on soils and vegetation. If a soil becomes more acidic, 
leaching of nutrients in the top soil layer (humus) is accelerated (ref. 7) , 
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This leaching affects, among other things, the rate of growth of trees; 
for example, forest productivity in Scandinavia has decreased bv one percent 
per year during the last few decades with high incidents of acidic preci- 
pitation (ref, 9). Ammonia (NH 3 ) plays an important role in neutralising 
acidic sulfur species. Consequently, if the global sulfur budget increases 
more rapidly tf in the global ammonia budget, a further increase in acid 
precipitation would be expected. 

Sulfurous Acid Radical CHSO 3 — mixing nitio 1.0 x 10"^ h 

It is of utmost importance to find out how gaseous sulfur is eventually 
converted into aerosol particulates. One reaction scheme that eventually 
produces H 2 SO 4 is via the IISO 3 radical. This HSO 3 to H 2 S 04 conversion 
is not immediate, however. It is therefore of considerable interest to 
discover the ultimate fate of the HSO 3 radical. Some complex reaction 
schemes for HSO 3 have been suggested (refs. 41-4,5), but none is supported 
by rate constants, 


Intermediate Sulfur Species 


Tl\e five intermediate sulfur species listed below are all very short lived 
and are therefore assumed to bo in photochemical equilibrium (PCB). Their 
lifetimes vary from 10 ““ to 10 ""* seconds and, hence, they react almost 
instantaneously to form other species. Due to their reactivity they cannot 
bo measured in situ. Their vertical profiles are calculated from known 
kinetic data. 


(i) 

Sul foxy 

(SO) ; 

(ii) 

Carbon sulfide 

(CS); 

fiii) 

Thiohydroxyl radieal(llS3 ; 

(iv) 

Sulfur trioxide 

CS03)5 

Cv) 

Atomic sulfur 

(SC^p)); 


tmixing ratio = 4.0 x 
(mixing ratio = 6.0 x 
(mixing ratio = 4.0 x 
(mixing ratio = 4,0 x 
(mixing ratio = -.0 x 


10 - 1 ^) 

10 " 20 ) 

10"‘l) 

10 ""“^) 

10 “^^) 
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ONB“DIMnNSIONAL TROPOSPHERIC PHOTOCHEMICAL MODEL 


The one-dimonsional, global, tz’oposphei'ic, photochemical model that 
has been under development at NASA/LaRC now contains a total of 108 chemical 
reactions and 12 photolytic reactions. Of this, the recent inlcusion of the 
sulfur group contributes 45 reactions (reaction 64 to reaction 108). Table 1 
lists the photodissociation processes that are included with their rates at 
the surface for a solar zenith angle of 45®. Table 2 lists all the chemical 
reactions with appropriate rate coefficients. 

Table 1, Photochemical reactions. 


Photodissociation 

No. Process Rate (Z°0) , s"^(45“) 


1 

O 3 + hv 0(^D) + Og 

6.3 E -6 

2 

O 3 + hv 0 + O 2 

3.0 R-4 

3 

NO 2 + hv ■+• NO + 0 

7.6 E-3 

4 

NO 3 + hv »>• NO 2 + 0 

8.2 E-2 

5 

NO 3 + hv ^ NO + O 2 

2.0 E-3 

6 

N 2 O 5 + hv NO 2 + NO 3 

2.0 E-5 

7 

HNO 3 + hv ->■ OH + NO 2 

1.6 E-7 

8 

HNOo + hv OH + NO 

4.8 E-4 

9 

H 2 O 2 hv -»■ 2 OH 

1.6 E -6 

10 

CH 2 O + hv HCO + H 

2,0 E-5 

11 

CH 2 O + 1\2 H 2 + CO 

6.5 E-5 

12 

CH 3 OOH + hv -)• CH 3 O + OH 

1.6 E -6 
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Table 2, Chemical reactions. 


Reaction 

No. 

Reaction 

Rate nxprossion 
(mole — cgs units) 

Reference 

1 

CH 4 + OH CH 3 + HgO 

2.4B-12*exp(-1710.0/T) 

Ca) 


CH 4 + 0 (^D) CH 3 + OH 

1.3E-10 

Ca) 

3 

dig + O 2 + M **" CH 3 O 2 + M 

2.20-31*(300.0/T)2'2 

Ca) 

4 

CH 3 O 2 + HO 2 CH 3 OOH + O 2 

3.0E-ll*oxpC“500,0/T) 

Cb) 

5 

CH 3 OOH + OH CH 3 O 2 + H 2 O 

6.2E-12*oxpC-750.0/T) 

Ca) 

6 

CH 3 O 2 + NO “*• CH 2 O + HNO 2 

0,0 

Cc) 

7 

CH 3 O 2 + NO 2 -»• CHoO + UNO 3 

0.0 

Cd) 

8 

CH 3 O 2 + NO CH 3 O + NO 2 

8.0E-12 

Ca) 

9 

CH 3 O + O 2 CH 2 O + HO 2 

S.0E-13*oxp(-2000.0/T) 

Ca) 

in 

CII 2 O + OH -V HCO + H 2 O 

1.7E-ll*pxpC-100.0/T) 

Ca) 

11 

HCO + O 2 CO + hOa 

5.0E-12 

Ca) 

12 

CO + OH -»• CO 2 + H 

1.35E-13ni - 

Ca) 

13 

CO + HO 2 CO 2 + OH 

l.OE-20 

Ce) 

14 

1I2 + OH ■+ H + H 2 O 

1.2E-ll*expC-2200.0/T) 

Ca) 

IS 

H + O 2 + M HO 2 + M 

5.5E-32*(300.0/T)l-‘* 

Ca) 

16 

0(^0) + H 2 O 2 OH 

2.3E-10 

Ca) 

17 

0 (^D) + H 2 H + OH 

9.9E-11 

Ca) 

18 

HO 2 + OH “>■ H 2 O + O 2 

4.0E-11 

Ca) 

19 

OH + Oil -> H 2 O + 0 

1.0E-ll*expC-S00.0/T) 

Ca) 

20 

HO 2 + HO 2 N 2 O 2 + O 2 

2.SE-12 

Ca) 

21 

H 2 O 2 + OH HO 2 + H 2 O 

1.0E-ll*expC-750.0/T) 

Ca) 

22 

IIO 2 + NO OH + NO 2 

5.4E-12*exp(2S0.0/T) 

Ca) 

23 

HO 2 + NO 2 HNOs + O 2 

3.0E-14 

Cf) 


10 


Table 2. tcontinued) . 


Reaction 

No. 

Reaction 

Rato E.xpression 
Cmole — cgs units) 

Reference 

24 

OH + NO 2 ^ HNO 3 

2.6B-30*(300/T)^-^ 

(a) 

25 

OH + NO ^ HNO 2 

2. OB-12 

(g) 

26 

HNO 3 + OH NO 3 + H 2 O 

S.5E-14 

(a) 

27 

HNO 3 rainout 

l.OB -6 

Cd) 

28 

NO 2 + O 3 NO 3 + O 2 

1.2E-13*exp(2450.0/T) 

(a) 

29 

NO + O 3 ">■ .NO 2 + O 2 

2.3E-12*exp(-l450.0/T) 

Ca) 

30 

,N 03 + NO 2 ^ N 2 O 5 

3.8B-12 

(g) 

51 

NO 3 + NO 2 NO + NO 2 + O 2 

2.3E-13*expC-1000.0/T) 

00 

32 

NO 3 + NO 2 NO 2 

8.7E-12 

(i) 

33 

NO + NO 2 + H 2 O 2 HNO 2 

6.0E-37 

(j) 

34 

N 2 O 5 ^ H 2 O -»■ 2 MNO 3 

l.OF-20 

(k) 

35 

N 2 O 5 NO 3 + NO 2 

5.7E-14*expC-10600.0/T) 

Cg) 

36 

0 + O 2 + M O 3 + M 

6.2E-34*C300.0/T)2- 1 

(a) 

37 

0(^D) + M 0 + M 

3.2E-11 

( 1 ) 

38 

OH + 03 "^ HO 2 + O 2 

1.6E-12*expC-940.0/T) 

Ca) 

39 

HO 2 + O 3 •> OH + 2 O 2 

l.lE-14*expC-580.0/T) 

Ca) 

40 

OH + HNO 2 -»■ NO 2 + H 2 O 

2.1E-12 

Cm) 

41 

OC^D] + CII 4 CH 2 O + H 2 

1.4E-11 

Ca) 

42 

N 2 H 4 + H N 2 H 3 + H 2 

9.9E-12*expC-1200.0/T) 

00 

43 

N 2 H 3 + N 2 H 3 -»■ 2 NH 3 + N 2 

K 43 << K 44 

Co) 

44 

N 2 H 3 + N 2 H 3 N 2 II 4 + N 2 H 2 

6.0E-11 

Co) 

45 

N 2 H 3 + H -> 2 NH 2 

2.7E-12 

Co) 

46 

1 NH 3 + 0 NH 2 + OH 

6.6E-12*exp(-3300.0/T) 

00 
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Table 2. (continued), 


Reaction 

No. 

Reaction 

Rate Gxpi’ession 
(mole — cgs units) 

Reference 

47 

NHa + 0(^0) NHa + OH 

2.5E-10 

(a) 

48 

Nils + OH NHa + HaO 

2.3E-l2*e.xpC-800.0/T) 

(n) 

49 

NII 3 + H *>• NHo + Ha 

l.OE-16 

(n) 

50 

NHa + 0 HNO + H 

1.8E-12 

fn) 

51 

NHa + 0 NH + OH 

1.8E-12 

(n) 

52 

NHa OH NII 3 + 0 

l.OE-13 

(n) 

53 

NHa + Ha NH 3 + H 

l.OE-16 

(n) 

54 

NHa + HO X- Na + HaO 

2.1E-11 

Cn) 

55 

NHa + HHa X- N 2 H 4 

l.OE-lO 

Co) 

56 

NH + NO X- Na + 0 + H 

4.7E-11 

Cn) 

57 

NH + Oa X- NO + OH 

6.0E-13 

Cp) 

58 

NH + NO Na + OH 

3.9E-11 

Cq) 

59 

HNO + On X- NO HOa 

2.1E-20 

Cn) 

(n 

HNO + M X- NO + H + M 

5.0E-08*expC-24500.0/T) 

Cn) 

61 

HNO + H NO + Ha 

5.0F-14 

Cn) 

62 

HNO + HNO X- NaO + HaO 

4.0E-15 

Cn) 

63 

H + H + M X- + M 

8.3E-33 

Cn) 

64 

HaS + OH X- MS + HaO 

l.C5E-ll*expC-200.0/T) 

Ca) 

65 

HS + 0 SO + H 

1.6E-10 

Cn) 

66 

HS + O 2 -X SO + OH 

l.OE-13 

Cn) 

67 

HS + NO X- products 

l.OE-12 

Cr) 

68 

SO + On SOn + 0 

7.5E-13*exp(-3250.0/T) 

Cn) 

69 

SO + NOn X- SOa + NO 

1.5E-11 

Cs) 

70 

CSa + 0 -X SO + CS 

3.1E-ll*exp(-640.0/T ) (£70 

= 0 . 8 ) (a) 
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Table 2, (continued). 


Reaction 

.Mo. 

Reaction 

Rato nxpression 
(mole — cgs units) 

Reference 

n 

CS 2 + 0 ^ S + COS 

3 . IR- 1 1 *e.xp ( -U4 0 , 0/T) ( f? 1 = 

' O.U(a) 


CSg + 0 ">■ St + CO 

3.1E“ll*exp(-640.0/T)(f7o s 

= 0 . 1 ) (a) 

?s 

CS + 0 -V S + CO 

2.2E-11 

(n) 

74 

CS + Ot so + CO 

•^S.OE-IS 

(s) 

75 

CS + Ot cos + 0 

<3.0E-18 

(s) 

7b 

s + Ot so + 0 

2 

(t) 

77 

SOt + OHC+M) HS03(+M) 

(8.2E-13*M)/(7.9E17 + M) 

(u) 

"8 

SOt + IIO 2 -»• SO 3 + Oil 

2.0E-1? 

(V) 

79 

sot + CHsOr. SO 3 + CH 30 

<3.3E-15 

(w) 

SO 

SO 3 + HtO HtSOi^ 

9.1E-13 

(n) 

81 

SO 2 + 0 so + 02 

(2 . 1E-10/.J,0 . S) *exp (-9980 . 0/T) (n) 

82 

11 s + 11 S + Ho 

2.5E-11 

(n) 

83 

S + II 2 11 s + H 

2 , 2E-25 

(n) 

84 

HoS + H HS + Ho 

1.29E-ll*exp(860.0/T) 

(n) 

85 

cos + H HS + CO 

2.2E-14 

(n) 

Sb 

COS + OH HS + COo 

5.7E-14 

(a) 

87 

CSo + OH HS + COS 

1.9E-13 

(a) 

88 

S + CSo So + CS 

6.5E-13 

(n) 

89 

S 0 - COS So + CO 

2.8E-12*exp(-2050.0/T) 

(n) 

90 

so + SO "*■ S + SOo 

3. on- 15 

(n) 

91 

SO + SO 3 2 SOo 

2.0E-15 

(n) 

92 

SOo + CHsC+M) CHaSOat+M) 

3.0E-13 

(n) 

93 

HS + HS ->• HoS + S 

1.2E-11 

(n) 

94 

CH3SH + 0 CII3 + HSO 

1.9E-12 

(X) 
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Table 2. (concluUed). 


Reaction 

No. 

Reaction 

Rate n.xpression 
(mole ~ cgs units) 

Reference 

95 

CH 3 SII + 

0 CH 3 SO + H 

1.9B-12 

t.t) 

96 

CH 3 SM + 

0 + CII 3 SOII 

1.9E-12 

(X) 

97 

CII 3 SH + 

Oil ">• products 

3.4E-11 


98 

CII 3 SCH 3 

+ 0 CII 3 SO + CH 3 

4.8R-11 

( 2 ) 

99 

CH 3 SCH 3 

+ 0 CII 2 S + CH 3 O 

6.3E-11 

Cx) 

100 

CII 3 SCH 3 

+ on products 

6.08E-12*e.\p(134.0/T) 

(aa) 

101 

COS + 0 

SO + CO 

2.1E-ll’'e.xp(-2200.0/T) 

Ca) 

102 

II 2 O + 0 

HS + OH 

2.4E-I2*e.xp(-1300.0/T) 

Ca) 

103 

so + O 3 

SO 2 + O 2 

2.5E-12*o.xp(-l050.0/T) 

(n) 

104 

SQz + 0 

+ M *»■ SO 3 + M 

3.4E-32*exp(-ll50.0/T) 

(n) 

105 

HSO 3 + ( 

OH *>• SO 3 + II 2 O 

l.OE-’' 1 

(bb) 

106 

SO 2 washout 

3.8E.6*(10-Z)/10 

(bb) 

107 

I ISO 3 ' 

washout 

2.3E-5*(10-Z)/10 

(bb) 

108 

Il 2 SOi, -*■ 

washout 

2.3E-5*(10-2)/10 

(bb) 
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Depending on the lifetime of a species, it can either be considered 
in photochemical equilibrium (PCR) or it will have to undergo transport, 

The rapidly reacting species have their vertical profiles determined by 
chemistry alone, while the vertical profiles of the more long-lived species 
are determined by the combined effects of chemistry and eddy ti’ansport, 
utilising a time-independent (sometimes called steady-state) continuity 
equation, The details of the continuity equation are described in Appondi.x A. 
Numerically, the continuity equation is solved using a finite-difference 
form known as the central difference scheme (see Appendix B for a discussion 
of this) , The boundary conditions for the continuity equation can be 
either a specified number density or a specified flu.x, as detailed in 
Appendix C, With the recent inclusion of the sulfur gi'oup, the model 
now contains a total of 37 chemical species. Of this amount, 17 are 
calculated using the continuity equation: nitric oxide (.NO), nitrogen 

dio.xide (NO2) , ozone (03), nitric acid (UNO 3), methylhydroporoxy CCH3OOH) , 
ammonia (NH3), hydrazine CN2H4), two hydrazine derivatives CN2II3) and CN2H2'' , 
hydrogen sulfide ClIjS), sulfur dioxide (SO2), carbon disulfide (CSo) . carbonyl 
sulfide (COS), dimethyl sulfide (CH3SCH33, methanethiol (CH3SH3, sulfuric 
acid (H2SO4) and sulfurous acid radical (HSO3), For the rapidly reacting 
species the condition of PCG is imposed. This includes the following 
species: the hydroxyl radical (OH), hydroperoxyl radical (HOo), nitrogen 

trioxide (NO3), dinitrogen pentroxide IN2O5), nitrous acid (NHO2), the 
amino radicals (NH2 and ,NH) , methylperoxy radical (CH3O2), metho,\y radical 
(CH3O), formaldehyde (HaCO) , formyl radical (HCO) , atomic oxygen ( 0 (^p)), 
the excited oxygen atom ( 0 (^D)), atomic hydrogen (H), nitroxyl radical (HNO) , 
thiohydroxyl radical (HS), sulfoxy radical (SO), sulfur trioxide (SO3) , atomic 
sulfur (S(^p)) and carbon sulfide (CS) . Certain very long-lived and 
well-mixed atmospheric species are specified as input parameters. This 
group includes molecular oxygen (O2) , molecular hydrogen (Ha), carbon 
dioxide (COa) , carbon monoxide (CO), methane CCH4), water vapor (H2O) , 
and hydrogen peroxide (H2O2) . Nine species are photodisassociated 
in 12 different paths in the troposphere in addition to participating 
in chemical reactions: they are ozone (O3) (2 photolysis paths), 

nitrogen dioxide (NOa) , nitrogen trioxide (NO3) (2 photolysis paths), 
dinitrogen pentoxide (N2O5) , nitric acid CHNO3), nitrous acid 
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tHNOi), hydrogen peroxide methylhydroporoxy tClI^OOin and formaldehyde 

CCH 2 O) (2 photolysis paths'), The details of the photodissociation calcula- 
tions are given in Appendix D. In order to ensure internal self-consistency, 
a convergence criterion is imposed. This criterion is mot when there is 
agreement between two successive runs up to the sixth decimal place. A 
further description of the convergence condition follows in Appendix B, 

Most atmospheric chemistry studios are perfoi’med vdth one-dimensional 
(vertican models, whereas two- and three-dimensional models are usually 
used for studies of atmospheric flow, dynamics and circulation. Since 
atmospheric chemistry is initiated by photodissociation of various mole- 
cules, the vertical coordinate is the important one in atmospheric chemistry 
studies. One-dimensional models are known as globally averaged models, 
because all parameters in these models should ideally be globally averaged. 
Two-dimensional models, or tonally averaged models, employ an averaging 
technique whereby parameters are averaged according to latitudinal bands, 
usually on the order of 5° to 10°, These models are mostly used to study 
certain flow phenomena of tonal character. Two-dimensional models are 
usually formulated fi*om a phenomenological point of view rather than from 
first principles. The eddy mixing is usually developed using mixing length 
theot’y. Hence, despite increasing the model by one more dimension, the 
transport coupling is still unsatisfactory since it relies on the choice 
of eddy coefficients. Three-dimensional models have been under development 
during tl^e last decade but have yet to be used very extensively due to 
their nearly prohibitive cost. The three-dimensional models in existence 
today have mostly been used for dynamic studies. The most sophisticated 
three-dimensional model in existence uses a more nine chemical reactions 
and seven chemical species. It has been calculated that, for each time 
seven new species are included, the computational cost rises an order of 
magnitude. Also, the largest digital computers available still do not 
have adequate storage to fully describe vertical transport. Hence, three- 
dimensional models must still rely on parameterized vertical eddy diffusi- 
vities, a fact that is not fully appreciated even in the scientific 
community. In one-dimensional models, a lot of storage is available to 
incorpoi'ate vast amounts of chemistry, which have most of their interesting 
features in the vertical dimension. Furthermore, most chemical interactions 
are more easily studied with one-dimensional models. 
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ROvSUI^TS AND DISCUSSION' 


Introduction 

In this section the sources and sinks of the sulfur compounds ai’o 
examined. Some sulfur species are produced chemically as well as emitted 
anthropogenical ly. Loss mechanisms include chemical reactions and 
heterogeneous losses, i.e., rainout, wa.'^hout, and dry deposition. The 
chemical production and loss terms for each species listed below are 
arranged in order of importance. 

Sulfur Dioxide (SO2) 

Chemical production .— Tlie chemical production of SOj always involves 
SO (reactions 68, 69, 103, 90, and 911. Of those reactions, only the 
recombination of sulfo.xy with molecular oxygen, Og, (reaction 68) and 
sulfoxy with nitrogen dioxide, NO^, (reaction 69) ni*e of major importance. 

SO reacting with otone, O3, (reaction 105) is of minor importance, while 
the reactions of SO with itself (reaction 90) and* with sulfur trioxide, 

SO3, (reaction 911 are negligible. 

Chemical loss .— The most important gas phase loss is duo to the 
I'oaction with the hydroxyl r.adical, OM, (reaction 77), A minor loss 
mechanism is due to reaction with the methylpero.xy radical, CH3O2 
(reaction 79). Reactions with the hydi'operoxy i*adical, HO2, (reaction 78); 
atomic o.xygen, 0, in the presence of a third body (reaction 104); the methyl 
radical, CH3, (reaction 92); and with atomic oxvgon, 0, (reaction 81) are 
all negligible. The lifetime of SO2 based solely on gas phase chemistry 
is about 2.5 days. 

A nthropogenic emission .— SO2 is knoMi as a product of fuel combustion 
(ref. 54). It has also been detected in the exhaust of turbines (ref. 44) 
and of diesel engines (ref. 36) . The application of certain t>'pos of 
fex‘tili:ers also emits SO2 (ref. 10). 

Natural emission .—’ Sulfur dioxide has been found in the smoko fi’om 
erupting volcanoes (refs. .“2, 38). 

Heterogeneous loss .— The heterogeneous loss of SO2 has been estimated 
to be about 50 percent of the total loss, and is modeled similarly to that 



by Turco (ref. 5), with tho c.^i^eption that the NASA/LaliC model has the 
tropopauso at 10 km while the model of Turco uses 13 km for tropopause 
altitude. Inclusion of tho heterogeneous loss of S02 lowers its lifetime 
to slightly less than two days, At the lower boundary, a flux 
of 1.2'’S X 10^® molecules cm"“ s"^ was used together with a depositional 
velocity of 1.0 cm s**^. This results in a surface mixing ratio of 0.5 ppb, 

Tho resulting vertical profile is sho\m in figure 2. Tho solid line 
represents the standard profile while the dashed lino is for the case without 
heterogeneous loss. In both cases a soro flux condition was imposed at tho 
upper boundary. Vortical profiles have been measured by dost (ref. 451 
and by Georgii (ref. 22 ), In general, .hose measurements show a decreasing 
mixing ratio with altitude similar to tho calculations with the NA,5A/ I.aRC 
model. One notable exception is a profile by Georgii (ref. 22) over the 
Atlantic Ocean (dotted lino). In this case there is an increase in mixing 
ratio with altitude. This is consistent with tho tlieory that the oceans 
act as a sink for SO 2 (ref. 20) . Also incluaed are some measurements by 
Maroulis (ref. 34) and by Jaeschke et al» (ref. 4b), Tho surface mixing 
ratio of the vertical pi*ofilos by Geotg^* (ref. 22) and dost (ref. 45) is 
about an order of magnitude higher than those calculated by the NASA/kaRC 
model. This is not very sui'prising, since Georgii's and dost's measurements 
were made over western Hurope, an area which is known to luavo a high back- 
ground concentration of SO 2 . Depositional velocities at the surface have 
a range of 0,1 to 2.5 cm s"^ (ref. 47). A value of 1.0 cm s’^ has been suggested 
as an average depositional velocity (ref. 48). The resulting vortical profiles 
due to different values of the depositional velocity are shown in figure 3, With 
a sl^w depositional velocity, a very rapid increase in mixing ratio with altitude 
results, and conversely with a high depositional velocity the mixing ratio 
increases slowly with altitude. In all calculations, the downward flux 
of SOa was kept constant at 1,275 x 10^® molecules cm"^ s”^. There are 
numerous estimates of global SOg emission. For example, Katt (ref. 49) 
estimated that 77 Mt were emitted. Subsequent estimates by Robinson 
and Robbins (ref. 31) indicate that the annual global SO 2 emission is 
14b Mt. Of this amount, 70 percent comes from coal combustion, lb 
percent from petroleum fuel combustion, and the rest from various processes 
such as petroleum refining and smelting. Approximately 95 percent of all 
tho anthropogenically emitted sulfur compounds are in the foi'm of SO 2 (ref. 12). 
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A recent estimate of global land sulfur dioxide emission was made by 
Shinn and Lynn (ref. 10). They estimate anthropogenic emissions to be 
183 Mt; 33 Mt of that amount is in the form of fertilizers. Natux’al 
enimissions Chiogenic and volcanic) account for 21 Mt of sulfur dioxide. 
Hence, man's activities account for nearly 50 percent of all sulfur 
dioxide that is emitted into the atmosphere. Since the bulk of this 
figure involves coal combustion, a switch to greater use of coal as 
recently suggested is likely to rapidly increase man's contribution 
to t!v) total sulfur budget. 


Carbonyl Sulfide (COS) 

Chemical production . The chemical sources of COS are the reactions 
of carbon disul'^'ide, CSg, with the hydroxyl radical, OH, (reaction 87), 
the reaction of molecular oxygen, 0 , with carbon sulfide, CS, (x’eaction 72) 
and the reaction of molecular oxygen, 0, with carbon disulfide (reaction 71). 

Only reaction 87 is of major importance. Reactions 72 and 71 are only minor 
pathways for production of COS. 

Chemical loss . — The loss mechanisms of carbonyl sulfide are totally 
dominated by the reaction of COS to OH (reaction 86). All other chemical losses 
of COS, i.e., reaction with atomic oxygen (reaction 101), reaction with atomic 
hydrogen, 11, (reaction 85) and reaction with atomic sulfur, S(^p), (reaction 89) 
can be neglected. The calculated lifetime of COS based on the loss mechanisms 
mentioned above is about four months. This value agrees reasonably well with a 
value of 0.6 year reported by Sze and Ko (ref. 4) and 200 days estimated by Logan 
et al. (ref. 6), 

Anthropogenic emission .— COS, which is a very stable molecule, is formed 
by manmade processes, such as manufacturing of petroleum (ref. 33) and of synthetic 
fiber (ref. 50) . 

Natural emission .— Some of the natural processes emitting COS include 
erupting volcanoes (ref. 32) and forest fires (ref. 51). 

In the model calculations, a surface boundary value of 0.5 ppb was used 
together with a condition of zero flux at the tropopause. The resulting vertical 
profile is show in figure 4. The measurements of COS have all been at the surface 
and indicate that carbonyl sulfide is very well mixed. Sandalls and Penkett (ref. 23) 
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Figure 4. Calculated vertical profile of carbonyl sulfide (COS) 


obtained a value of 0,51 ppb, Torres et al. (ref. 25) measured an average value 

of 0.512 ppb over a wide range of latitudes, while Maroulis et al. Cref. 24) 

found an average value of 0.467 ppb over 3 different locations in the United 
States. Some early data by Hanst et al, (ref. 26) indicate a value of 0.2 ppb 

at the surface, but this measurement was associated with lar'ge e.xperimental 

uncertainties. The vertical profile calculated by the model shows a nearly 
constant mi,xing ratio in the troposphere as would be expected from its long 
lifetime, COS diffuses into the stratosphere, where it undergoes photolysis 
forming SO (ref, 1). Hence, COS would be an important precursor to the 
formation of the particulates that constitute the aerosol layer. 

Hydrogen Sulfide (HaS) 

Chemical production .-— The only knoim chemical production of IIjS in the 
atmosphere is due to the reaction of the thiohydroxyl radical, HS, with 
itself (reaction 93). This is only of very minor importance. 

Chemical loss .— The primary loss mechanism for hydrogen sulfide is the 
reaction of HaS to OH (reaction 64). The losses due to reactions with atomic oxygen 
(reaction 102) and with atomic hydrogen (reaction 84) are negligible. The calculated 
lifetime of HaS based on these three chemical losses is a little more than one day. 

Anthropogenic emission .— Some of the industrial activities that have been 
identified as emitting HaS include wood pulping (ref. 29) and sewage treatment 
(ref. 30). 

Natural emission .— HaS is formed primarily by microbial activity 
(ref. 27). Volcanoes also emit hydrogen sulfide (ref. 28) as does animal 
waste (ref. 8). 

The only measurements available for HaS are those of Geory.ii (ref. 22) 
and of Slatt et al. (ref. 52). The measurements by Georgii indicate a 
surface mixing ratio of 1.5 to 3.0 ppb. This is close to the estimate of 
0.2 ppb by Robinson and Robbins (ref, 31). Using this value as the lower 
boundary condition and zero flux at the tropopause, the model calculated 
the profile shown in figure 5. The model shows a very rapid decrease in the 
mixing ratio with altitude similar to Georgii's measurements. 
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Figure 5. Calculated vertical profile of hydrogen sulfide (H 2 S). 



Carbon Disulfide CCS^) 

Chemical production .-— No kno\vn chemical reaction in the atmosphere 
produces carbon disulfide. 

Chemical loss . — The primai'y loss mechanism is the reaction with the 
hydro.yyl radical, OH, (reaction 87), Minor loss contributions come from 
reactions with atomic oxygen (reactions 70, 71, and 72). The reaction of 
CS 2 to atomic sulfur, S(^p), (reaction 88 ) is negligible. The lifetime 
based on these reactions is slightly more than one month. This value compares 
reasonably well with the value of 0.2 years given by Sse and Ko (ref. 4). 

Anthropogenic emission .— CS^ is kno»vn to result from the manufacturing 
of petroleum (ref. 33) and of synthetic fiber (ref. 36). 

Natural emission . — Carbon disulfide has been detected in the smoke 
from volcanoes (ref. 32). Also, some measurements in seawater (ref. 37) 
seem to indicate a natural source of CS^ , 

Measurements of CS^ are very sparse and confined to one location in 
England (ref, 23). The average surface value obtained was 0.19 ppb. No 
vertical measurements are available. Using the value of 0.19 ppb as the 
lower boundary condition and zero flux at the tropopause, the model calculated 
the profile sho\m in figure 6 , CS^ decreases relatively slowly with altitude, 
as would be expected from its moderately long lifetime. 

Sulfuric Acid (H SO ) 

^2 it 

Chemical production . — Sulfuric acid is formed chemically by the reaction 
of sulfur trioxide, SO 3 , with water vapor, HgO, (reaction 80). 

Chemical loss .— No gas phase chemical loss mechanism has yet been 

identified for H„SO, . 

2 4 

Anthropogenic emission . Sulfuric acid has been detected in the exhaust 
of automobiles (ref, 36) and in the manufacturing pi*ocess of H^S0^^(ref. 39). 

Natural emission .— H^SO^^ is present in volcanic emission (ref. 38), 

Heterogeneous loss . — The only knom loss mechanism for sulfuric acid 
is heterogeneous loss. In the NASA/LaRC model we used a heterogeneous loss 
rate similar to that of Turco et al. (ref. 5), again with the exception of a lower 
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vertical profile of carbon disulfide (CS2}. 




tropopausc altitude (10 km vs. 13 km). The lifetime of H2SO4 at the 
surface based on this heterogeneous loss rate is about 0,5 day. 

It is of utmost importance to learn the fate of the H2SO4 molecules, 

It can occur in liquid form and rainout in the fox’m of acid ■ xin which 
affects the environment adversely, or it can occur in solid form, which 
affects the growth of the aerosol layer and, hence, ultimately will 
affect the climate. In the model, a tei’o flu.\* condition was imposed at 
both boundaries, The resulting vertical px’ofile is shown in figure 7 . 

There are no tropospheric measurements of H^SOi, for comparison. 

Sulfux’ous Acid Radical (USOa) 

Chemical production .— The only know chemical production is due to 
the reaction of SO2 with OH (reaction 77 ) . 

Chemical loss .— Gas phase loss occurs as a result of the reaction 
of HSO3 with OH (reaction 105 ), The homogeneous lifetime of HSO3 is 
appro, \imately 15 hours at the surface. 

Anthropogenic and natural emission . — No anthropogenic or natural 
emissions are presently know. 

Heterogeneous loss . The heterogeneous loss tei-m is modeled similar 
to HgSOi, (ref. 5 ). The combined lifetime due to both homogeneous and 
heterogeneous chemistry is about six hours at the surface. 

It is of great interest to find out the details of the fate of the 
HSO3 radical because it offers a loss mechanism for sulfur before oxidation 
to sulfuric acid occurs. As in the case of H2SO4 , no atmospheric measure- 
ments of HSO3 exist. The free acid H2SO3 is not believed to exist, wliile 
the bisulfates containing HSO3 are well know (ref. 53 ) , The boundary 
conditions chosen for HSO3 are similar to those for H2SO4. The resulting 
profile is show in figure 7 . Tlie profiles of HSO3 and H2SO4 follow each 
other in shape very closely. This is due to the fact that the dominating 
heterogeneous loss term is identical for both molecules, 



Intermediate Sulfur Species 

Introduction . -- The five short-lived intermediate sulfur species have 
only chemical production and loss mechanisms with the exception of sulfur 
trioxide, SO 3 . Anthropogenic emissions of SO 3 have been detected in fuel- 
combustion processes (ref. S4) and in the manufacturing process of sulfuric 
acid (ref. 5S). The following subsections will discuss the chemical 
production and loss mechanisms for each of the five intermediate species. 

Again, the chemcial reactions will bo mentioned in order of importance. 

Sulfo.xy C50) .— SO is formed primarily by oxidation of the thiohydro.'cyl 
radical, HS, with molecular oxygen, O 2 (reaction 66 ). Of minor importance 
are the oxidations of carbon disulfide with 0 (reaction 70), carbon sulfide, 

CS, with 0 (reaction 74), atomic sulfur [S(^p)j with O 2 (reaction 76) and 
carbonyl sulfide, COS, with 0 (reaction 101). The productions of SO due 
to the oxidations of the thiyl radical, US, with 0 (reaction 65) and of 
sulfur dioxide, SOg, with 0 (reaction 81) arc negligible compared to the 
other reactions producing SO. Losses are due primarily to the oxidation 
of SO with Og (reaction 68 ) . Minor losses are caused by the reaction of 
SO with NO 2 (reaction 69) and by SO reacting with ozone, Oj, (reaction 103). 
Negligible loss mechanisms include the reaction of SO with itself (reaction 90) 
and of SO reacting with SO 3 (reaction 91). Hence, the five loss mechanisms 
for SO nearly always produce SOg. In only one instance, the reaction of SO 
with itself (reaction 90), is SO 2 not produced. Thus, SO acts as an important 
precursor to SO 2 . The calculated lifetime of SO is on the order of lO"^ s. 

The vertical profile is plotted in figure 8 . Due to its short lifetime and 
high reactivity, measurements cannot be made. This is true for the four 
remaining intermediate sulfur species also. 

Carbon sulfide (CS) .— CS is formed primarily when carbon disulfide 
is oxidized by atomic oxygen (reaction 70) and, of negligible importance, 
when atomic sulfur reacts with CSg (reaction 88 ) . The two primary loss 
mechanisms are oxidation of CS with Og (reactions 74 and 75) . The reaction 
of CS with Og takes two paths with equal preference. A negligible loss 
for CS is the oxidation with 0 (reaction 73). The lifetime of CS based on 
these reactions is on the order of 10"^ s, well justifying the PCE assump- 
tion. The vertical profile is shorn in figure 8 . 
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Thiohydroxyl radical CHS) .«»° The thiohydi*oxyl radical is a highly roactivo 

species formed by seven r 'actions and destroyed by five . Three reactions are of 
importance for its formation: hydi’ogon sulfide I’eacting with the hydroxyl 

radical, OH, (reaction 64), carbon disulfide reacting with OH (reaction S**), 
and carbonyl sulfide reacting with OH (reaction sfO . Minor nroduction of 
HS is due to the oxidation of HaS with 0 (reaction 102} and the reaction 
of H with H 2 S (reaction 84), Of negligible importance are the reactions 
of COS with H (reaction 85) and S reacting with H 2 (reaction 85). Losses 
are due almost entirely to the oxidation of O 2 (reaction 66 ). Of negligible 
importance are the remaining four loss mechanisms: HS reacting with NO 

(a’eaction 67), HS oxidized by 0 (reaction 65), HS reacting with H (reaction 82), 
and HS reacting with itself (reaction 93). The resulting lifetime of the 
thiohydroxyl radical is very short, on the order of 10 "^ s. Its calculated 

vertical profile is plotted in figure 8 . 

Sulfur trioxide (SO 3 ) .— • As mentioned previously, sulfur trioxide 
is the only intermediate sulfur species that has nonchemical production 
terms, The relative importance of chemical vs. nonchemical contributions 

I 

to the total SO 3 budget is poorly understood. Cliemically, SO 3 is produced 
primarily by the reaction of HSO 3 with OH (reaction 105) and by SOg 
reacting with the methylperoxy radical, CH 3 O 2 , (reaction 79), Minor 
contributors to SO 2 production are SO 2 reacting with HQ 2 (reaction 78) 
and SO 2 being oxidized by 0 in the presence of a third body (reaction 104). 

The main path for SO 3 destruction is the reaction with water vapor, H^O, 

(reaction 80). Of negligible importance is the reaction of SO 3 with SO 
(reaction 91). Sulfur trioxide is very short lived. Its lifetime is on 
the order of 10"® s. The vertical profile due to the chemistry described 
above is plotted in figure 8 . 

Atomic sulfur [S(^p)] .— The production of atomic sulfur is due 
almost entirely to the reaction of CSg with 0 (reaction 71). Of negligible 
inportance are the four reactions: CS with 0 (reaction 73), SO with itself 
(reaction 90), HS reacting with H (reaction 83), and HS reacting with itself 
(reaction 93). Loss of S is primarily by oxidation with molecular oxygen 
(reaction 76). Minor contributors to loss of S are reactions with CS 2 
(reaction 88 ) and with COS (reaction 89). The reaction of S with H 2 
(reaction 83) is virtually negligible as a loss for atomic sulfur. The 
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lifetime of atomic sulfur is extremely shorti about 10”*' s. The vertical . 
profile of S is shouTi in figure 8. (Observe that the horizontal scale is 
shifted for S.) 

FUTURE PERTURBATIONS TO THE SULFUR BUDGET 

A couple of decades ago, the only global scale implication due to 
sulfur species was thought to bo the formation of stratospheric particulates 
(ref. 2). This argument was based on the relatively moderate lifetimes of 
most sulfur species, Lifetimes on the order of a few days to a weoK were 
thought to preclude any real important contributions to global ti’opospheric 
chemistry. We now know that sulfur species play an integi'al part in tlie 
highly coupled and complex chemical system that constitutes the troposphere. 
SO 2 emissions had an upward trend during the IDoO's. This trend was 
partially reversed during the 1970 's. Some actions responsible for this 
decrease were a greater use of clean-burning natural gas and of crude oil 
with low sulfur content. Also, a slower economic growth coupled with 
energy conservation contributed to a deci*ease in SO 2 omissions. Recently, 
this trend has started upward again, The global concentrations of sulfur 
compounds are vei*y likely to see a rapid increase in years to come duo to 
energy policies. Already emission standards for many oil-burning power 
plants and industrial plants are relaxed to allow for usage of lower grade 
crude oil containing more sulfur. Anthropogenic emissions are currently 
responsible for approximately half as much of the atmospheric sulfur 
concentrations as are the natural emissions, but some 20 to 25 years from 
now manmade emissions may equal those of nature and may actually at the 
turn of the century surpass nature's emissions in the northern hemisphere 
(ref. 12). The estimates of global emission of sulfur are widely varying. 
Katz (ref. 49) made some early estimates and calculated that 77 Mt 
.of SO 2 was emitted in 1943. The Study of Critical Envii'onmental Problems 
(ref. 56) established a global emission of 93 Mt/yr in 1967 to 1968 
Robinson and Robbins (ref. 31) estimated the annual total .SO 2 emission in 
the raid-1960' s to be 146 Mt. A more recent estimate (ref. 10) 
calculates the global manmade omissions of SOj to be 183 Mt. It is 
significant to note that of the total global emissions of SO?,, 93.5 percent 
occur in the northern hemisphere and only 6.5 percent in the southern 
hemisphere (ref. 56). The most recent estimate (ref. 10) of global 



bactoriogenic production of SOg is 210 Mt. Fertil iters are thought 
to account for 33 Mt of the manmade emissions. This portion of the 
sulfur budget is very likely to increase further. In the SCEP report 
(ref, 36), some predictions for the future wore made. An annual 
growth in fossile fuel usage of 4 percent between 1970 and 1980 was 
assumed, with a 3.S percent growth rate from 1980 to 2000. By the year 
2000 A.D.f a global anthropogenic omission of 2?S Mt/yr would occur 
(ref. 56). However, at the present time it seems that any prediction 
of future crude oil is risky at best, 


It is very evident that future sulfur levels in the atmosphere are 
a strong function of man's activities. Many questions concerning 
sulfur species remain to be answered: for example, the exact amounts 

of naturally vs. anthropogenically produced sulfur need to be ascertained. 
Also, what is the global concentration of SOa? Present measurements are too 
inconclusive and in part contradictory. Furthermore, short-term and long-term 
impacts of SOg omissions on the sulfur budget must bo elucidated, The primary 
chemical loss for SO2 is OH. The hydro.xyl radical (OH) is very reactive 
and acts somewhat like a "tropospheric vacuum cleaner." In figure 9, the 
resulting SOg profiles due to a doubling of the current OH level and a 
decrease by 50 percent in the present OH level fire plotted. Thus, any 
increase in the SO2 level would compete with many other species for the 
OH radical, The budgets of these competing trace gases need to be addressed. 
The importance of homogeneous vs, heterogeneous loss for SO2 needs further 
research. Presently, it is believed that 50 percent is lost homogeneously 
and SO percent is lost heterogeneously, but this breakdown is a very crude 
estimate, The depositional velocity of SO2 needs further pinpointing. 
Estimates vary from 0.1 to 2.5 cm s"^ with 1.0 cm s"^ as an average. The 
vertical profile of SO^ is very dependent on this quantity, as already 
pointed out in figure 3. It is also urgent to find the mechanism that 
controls the SO2 to sulfate (gas-to-particle) conversion. From figure 1 
we can see that SO2 is converted to H2SO4 via the intermediate species SO3. 

A possible alternate reaction path involves the HSO3 radical. The fate of 
this molecule is virtually unknom. Also, the kinetic data concerning the 
reactions that convert SO2 to SO§" need to be improved. Finally, the 
role of ammonia, NH3, in aei'osol formation needs to be addressed. NH3 is 
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largely responsible for neutralising the acidic sulfur species. Hence, if 
the sulfur budget increases more rapidly than the ammonia budget, an 
increase in ar-idic precipitation over the already high levels would result. 

In conclusion, we have seen that the sulfur family, because it is 
largely so anthropogenic plays an increasingly important role in tropospheric 
chemistry. Additional research is needed to answer many of the questions 
related to the tropospheric sulfur budget. The answers to tliese questions 
will have a direct impact on many fields of scientific study, e.g., 
radiative transfer, climatology, meteorology, agronomy, ecology, and 
limnology, just to mention a few. 
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APPENDIX A 


THE CONTINUITY EQUATION 


The continuity equation can bo written in either number density form 
or mixing ratio form. In the NASA/LaRC one-dimensional tropospheric 
photochemical model, the mixing ratio form is used. The vertical profile 
of a long-lived species is then e.xpressed as: 

9<j>. 


where is the vertical flux (molecules cm"^ s"*) of the ith species, 

Q. (n.) are the chemical production terms, and L.(n.)Mf. are the chemical 

1 J 1 J 1 

loss tei'ms of the ith species; M is the total number density (molecules cm“^), 
f^ is the mixing ratio of the ith species, and M and f are related 
by the e.xpression: 


where n^ is the number density of the species under consideration. 

The vertical flux of the ith species, can be written in terms of a 

parameterization, K„ : 


<j). = -K M 
1 z 


'_J 

3z 


(A-3) 


The term K, is an empirical constant (cm^ s“^3 usually called the 
eddy diffusion coefficient. The word "eddy" is somewhat of a misnomer 
since it indicates that the diffusive process occurs on a small scale. 

On the contraxy, most vertical transport in the troposphere takes place 
on a vexy large scale. Substitution of equation (A-5) into equation (A-1) 
yields : 


K. M 


af, 

Jz 


= -Q. (n . ) + L. (n . ) Mf . 

^ j 1 ^ j 1 


(A-4) 
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Because equation (A-4) usually depends on many species othei* than the ith, 
it is a highly coupled, nonlinear, partial differential equation that has 
to be solved numerically. 


In the case of short-lived species, the chemistry dominates the vertical 
distribution. For these species the vertical transport terms [i.e. the 
terms on the left-hand side of equation (A-4)] can be neglected. This 
condition is knoivn as photochemical equilibrium (PCE) . If the PCE 
assumption is justified, the solution of equation (A-4) is simplified con- 
siderably since we can solve e.xplicitly for f^: 


^i " L.(n.)M 


(A-5) 


Hence, in the cases where PCE can 
efficiency, both in terms of time 


be used, a high degree of computational 
and money, can be reached. 
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APPENDIX B 


NUMERICAL SOLUTION OF THE CONTINUITY EQUATION 

The NASA/LaRC one-dimensional, global, troposphei'ic, photochemical 
model calculates the continuity equation in the mixing ratio form, 
previously given in equation (A-4) : 




= -Q.(n.) + L.(n.) Mf. 
13 3 -^ 1 


By expanding on the derivative with respect to altitude and rearranging the 
terms in equation (B-1) we obtain: 

3^fi 9 

M — ^ CK, M) -i - L. (n.) Mf. + Q. (n.) = 0 

“ a - 2 a, a? •'•J J-J 


After dividing this equation by M we obtain: 


'i , 1 3 nr 


(K M) - L,(n,)f, + 


Q, CnJ 


1.1 _ 


M 8z '• z 3z “i^‘‘j^^i M 


Equation (B-3) defines a system of equations for each species, f^. 
A finite-difference scheme is used to solve the equations (ref. 57) 
The finite-difference equations are defined as: 


a2 f _ 4. £ 

o ^lv + 1 ^ ^k— 1 

^ ^ (for second order equations) 

2 A 2 


. ^k.i - ^k-i 


(for first ox'der equations) 
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Substituting these equations into equation (B-3), rearranging terms and 
dropping the species index i and subscript z on K, , we obtain: 


k+1 


’ Kk A ’ 



* ^k.i 

1 

< I 

1 

2As 

k 

Az^ 

[Az^ 2AzJ 


M 


(B-6) 


1 HK M) 

The term A represents: |;| — ^ , and the index k refers to the 

spatial derivative z. Equation (B-6) can be written in finite-difference 
form as: 


A"f +’5*f +^*f =D 
^ ^k+1 ^k ^ \-l 


(B-7) 
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APPBNDIX C 


BOUNDARY CONDITIONS 


Two types of boundary conditions are available to solve the continuity 
equation. Wo can either specify a mi.xing ratio f^ or a flu.x 
upward or downwai’d at the upper and lower boundaries. In general, the 
boundary condition can be written as: 




CC-1) 


If we specify a ini.xing ratio f^ at the lower boundary, a^ =0, b^ = 1, 
and equation (C-l) becomes: 


f. 


1 


c. 

X 


CC-2) 


or, in finite-difference foriiU 

h • = h 


CC-3) 


where the subscript ^ is used to denote lower boundary. If a mixing 
ratio is specified at the upper boundary, a. = 1, b. = 0, and equation 
(C-l) becomes 


a. 

1 



c 


i 


CC-43 


In finite-difference form (with subscript u for upper boundary) : 


B • f + ^ • f = 
u u u u u 


This result is obtained from an expansion of equation (C-4) 


(C-5) 
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f - f , 
u u-1 


» c. 


CC-0) 


This can be rearranged to: 


u 


u-1 




CC-7) 


which in finite-difference form can be written as equation (C-5) , with 

^ = 1, & =0, and 5 n c.. At. 
u u ’ u u 

For an upward or downward flux across the upper boundary, the boundary 
condition reduces to: 


3z 


= c(<t>3 


u 


The finite-difference form looks similar to equation (C-5) 


CC-83 


• fW„ + CC-9) 

In this case, = 1, = -1> *^*-‘**^u “ 

For a flux, upward or do\mward, across the lower boundary, equations CC-8) 
and (C-9) are still valid, except in this case = -1, = 1, 

= CC4>)j, Az, and A(<|i)£ = 0. 

Equation (B-7) from Appendix B, together with equations (C-3) and (C-53, 
forms a so-called "block tridiagonal" system of equations given by: 
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There are several numerical techniques available to solve the tridiagonal 
matrix described by equation (C-10). In the NASA/LaRC model, a Gaussian 
elimination method without pivoting is used (ref, 57). 
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APPENDIX D 


PHOTODISSOCIATION RATE CALCULATIONS 


In ordei‘ to calculate the photodissociation I'atos, one first must 
determine the amount of incoming solar radiation. The incident I'adiution 
is a function of wavelength X, altitude z, and solar zenith angle 
0 [i.e., I » ICXiZ|0)]. The expression for incident solar radiation 

is due to Leighton (ref, 58): 


I(X,z,03 a lQ(X)*exp[-To^CX)sec 6] *joxp[-T^^CX) - t^,(X)] sec 0 
+ fl - exp(-TjjCX) - T^(\)]sec ^‘cos 0 


(D-1) 


where I^t-X) is the incident solar flux at the top of the atmosphere. 

This has been tabulated by Ackerman (ref. 59). The term exp[-To (X) 

3 

• sec 0] is the attenuation due to ozone absorption. The next term 

exp[-Tj^(X) - T^,(X)]sec B is the direct solar attenuation due to aex’osol particle 

scattering (Xp) and Rayleigh scattering (t^) . Tlie last term [1 - e.\p[-r^^(X) 

- T^.(X) sec 0] cos 0 is the diffuse solar radiation attenuated by aei'osol 

and Rayleigh scattering. Values for t (X) and x (X) have been 

p r 

tabulated by Elterman (ref. OO) . Once the incident solar radiation, 

I(X,z,0) is known, the photodissociation rates, J^, can be calculated 
according to: 


J.(X,z,0) = EAI.CX,z,0)cr.(X) (D-2) 

^ X=1 ■' ^ 


where <^j^CX) is the molecular cross section of the ith species. The 
solar zenith angle, 6, is calculated as: 

cos 6 cos 'P * cos & * cos t + sin (}) * sin d, (D-.5) 

where i{i is the latitude, d is the solar declination angle and t 
is the local hour angle of the sun. For a given latitude (tj) and solar 
declination (d), the local hour angle of the sun (t) varies from -180° 
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to 180° 
noon is 


with each hour corresponding to a 15° increment, i. 
0°. 11 a.m. is -15°, 1 p.m. is +15°, etc. 


e., local 
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APPENDIX E 


CONVERGENCE CRITERIA 


Tlje model imposes two separate convergence criteria, one for the 
species in photochemical equilibrium and one for the species that are 
calculated using the continuity equation. The number of iterations 
needed is a strong function of the initial profiles that are prescribed, 

In general, the flow of the calculations of the model is as follows: 

First, the initial profiles for all species are prescribed. Secondly, 
the reaction rate constants are calculated. Next^ the incident solar 
radiation is calculated, and based on these results, the photodissociation 
rates are computed, The model then calculates the vertical profiles 
of tiie short-lived species in photochemical equilibrium. Finally, the 
long-lived species that are transported are calculated. The vertical 
profiles that are obtained for the species are compared to the previous 
iteration and recalculated until the convergence criteria is achieved. 

For the short-lived species, the convergence criteria is 


1 


Uii)2 


10' 6 


CE-n 


and for the transported species the criteria is 


1 


Cn.), 

1 1 


< lO-** 


(E-2-) 


Computationally, the photochemical equilibrium species are calculated 
much more rapidly than the transported species; therefore, a somewhat 
more stringent convergence criteria can be imposed on the PCE species. 
With reasonably close initial guesses of the vertical profiles (i.e., 
vertical profiles within a couple of orders of magnitude of the final 
profiles) , convergence is achieved after three to four iterations for 
the short-lived species and five to six iterations for the long-lived 
species. 
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